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Chromatin is a complex of DNA and spe- 
cific proteins forming an intermediary level 
of organization of eukaryotic genomes, be- 
tween double-stranded DNA and chromosome. 
Within a generic modeling of the chromatin as- 
sembly, we investigate the interplay between 
the mechanical properties of the chromatin 
fiber and its biological functions. A quantita- 
tive step is to relate the mechanics at the DNA 
level and the mechanics described at the chro- 
matin fiber level. It allows to calculate the com- 
plete set of chromatin elastic constants (twist 
and bend persistence lengths, stretch modulus 
and twist-stretch coupling constant), in terms 
of DNA elastic properties and geometric fea- 
tures of the fiber. These elastic constants are 
strongly sensitive to the local architecture of 
the fiber and we argue that this tunable elas- 
ticity might be a key feature in chromatin func- 
tions, for instance in the initiation and regula- 
tion of transcription. Moreover, this analysis 
provides a framework to interpret micromanip- 
ulations studies of chromatin fiber and suggests 
further experiments involving intercalators to 
scan the tunable elasticity of the fiber. 



Introduction 

Chromatin is a protein-DNA complex forming the 
basic material of chromosomes of all eukaryotic organ- 
isms and strikingly conserved during evolution jjj]. It 
provides an intermediary level of organization between 
the underlying DNA double-helix and the whole cell 
machinery. Chromatin achieves essentially two funda- 
mental functions: it compacts efficiently DNA into the 
ceil nucleus, and it compacts it operationally in or- 
der to allow targeted gene expression. We propose a 
mechanical approach of these two questions. It sug- 
gests that the chromatin structure might be selected 
according to its elastic properties, which might be ex- 
ploited during transcription by the various enzymes at 
work during this process, for instance to locally con- 
dense (gene silencing) or decondense (gene activation) 
the fiber. 



Chromatin assembly 

We introduce a basic geometric modeling of chro- 
matin assembly, which allows to incorporate acknowl- 
edged data on DNA structure, its bend and twist per- 
sistence lengths, and data on the nucleosome struc- 
ture into an analytically tractable framework. Chro- 
matin is composed of a double-stranded DNA molecule 
wrapped from place to place around protein cores (hi- 
stone octamers). Focusing on the native case where 
these cores are fixed on the DNA, we model separately 
the nucleosomes, i.e. the DNA-wrapped histone cores, 
and the linkers, i.e. the naked DNA segments connect- 
ing the nucleosomes. 

Since we are looking for generic properties of the 
chromatin fiber, originating from its assembiy, se- 
quence effects are ignored. We suppose that all the 
linkers have the same number n of base pairs; it corre- 
sponds to currently observed phased nucleosomes ||. 
A posteriori, we may argue that self-organization of the 
fiber architecture leads to such a regular positioning of 
the nucleosomes in order to get proper structural and 
elastic properties. In fact, we need this hypothesis to 
be satisfied only locally, over a few linkers. Linker DNA 
is described as a non extensible semi-flexible polymer 
within the continuous worm-like-chain model supple- 
mented with twist elastic degree of freedom |^] ; it thus 
involves two persistence lengths, respectively of bend 
(A — 53 nm) and twist (C = 75 nm), the given values 
corresponding to lOmM NaCl Linkers are straight 
in absence of applied constraints. 

Nucleosome structure is now well-known thanks to 
high resolution cristallography data J5|. The net effect 
of a nucleosome on the DNA trail can be described as a 
rigid kink connecting two successive linkers. The slope 
of the DNA path on a nucleosome being fixed (angle 
4.5°), a single angle $ between entering and outgo- 
ing linkers thoroughly prescribes the geometry. This 
picture allows an effective modeling of various physico- 
chemical effects, for instance electrostatic repulsion be- 
tween linkers and its variable screening when the ionic 
strength is increased. It accounts also for histone HI 
binding in the neighborhood of linker entry/exit points 
as well as histone tails influence and its variation upon 
posttraductional covalent modifications of these tails 
Iq @ l§ ||' We here ignore interactions between nu- 
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cleosomes, which is valid as soon as the internucleoso- 
mal distances are larger than the interaction range; in 
any cases, the model describes the geometric contribu- 
tion to chromatin elastic behavior originating from the 
DNA elastic properties. 

Our model, quite similar to the so-called two-angle 
model of Woodcock et al. Jh]], amounts to build the 
chromatin fiber nucleosome after nucleosome. The ro- 
tational positioning of a nucleosome with respect to 
the previous one is entirely prescribed by the twist an- 
gle of the DNA linker connecting them. In the relaxed 
state, it is equivalently prescribed by the linker length 
n (in bp) since the twist angle is then precisely equal 
to r = 2n n/10.6 (recall that there is 10.6 bp per DNA 
turn) . Each step of the assembly is described analyti- 
cally, but we implemented this model within a Maple 
program in order to handle an arbitrary number of nu- 
cleosomes and to perform a quantitative analysis of the 
chromatin fiber geometric properties. 




Figure 1: Varying the linker length n and the en- 
try/exit angle $ recovers all the structural diversity 
observed experimentally, for instance ribbon ($ = 90°, 
n = 58 bp), columnar structure ($ = 50°, n = 42 bp) or 
double-helical structure ($ = 60°, n = 56 bp). 

Structural properties 

The implementation of the above model straightfor- 
wardly gives all the geometric properties of the chro- 
matin fiber. Homogeneity of the local architecture 
(uniform values for n and $) enforces a helical symmet- 
ric structure for the fiber; we shall speak of "chromatin 



superhelix" (SH), and denote A its axis. A complete 
study of the symmetry properties and quantities rele- 
vant to the mechanics of the fiber has been presented in 
a previous paper jn|. We here only underline the high 
structural sensitivity of the fiber, evidenced on Figure 
1: for various values of n and we obtain the whole 
variety of structures suggested from experimental stud- 
ies m Q Jli). These qualitative results support 
our minimal two-parameter modeling, since it is suffi- 
cient to recover the actual diversity of chromatin struc- 
tures. A noticeable point, evidenced on Figure 2, is 
the fact that the diameter remains about 30 nm what- 
ever the structure; the so-called 30nm-fiber observed 
experimentally might thus correspond to far different 
structures. Excluded volume has to be taken into ac- 
count to discard unrealistic structures; we determined 
a sufficient condition to avoid steric hindrance, which 
is shown on Figure 2. 

One of the conclusions of our structural study is the 
fact that the connection between the microscopic pa- 
rameters and the SH geometric characteristics is too 
complex and multivariate to get small-scale informa- 
tions from structural observations of the fiber; we now 
turn to the analysis of its mechanical properties. 




linker length (bp) 



Figure 2: Diameter of the chromatin fiber (upper blue 
curve) and its pitch V (lower red curve). The thin 
green line gives a lower bound V c for the pitch above 
which no steric hidrance occurs within the fiber, in 
constrast to the dashed regions where V < V c . 

Elastic properties 

In order to obtain the leading order of the elastic 
constants of the fiber, we merely consider a straight 
fiber with cylindrical symmetry. Describing the elastic 
properties of the chromatin fiber is thus nothing but a 
problem of spring mechanics [^. Nevertheless, the ar- 
chitecture of this "spring" is much more complex than 
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a simple helical coiling and we expect that the detailed 
structural features of the chromatin assembly strongly 
influence the behavior at the fiber scale. Our modeling 
allows us to compute analytically the elastic constants 
describing the linear response of the fiber when a force 
F and a torque M are applied at its ends, hence to 
investigate quantitatively how the elastic response of 
the fiber varies with its relaxed structure. Note that 
this structure is itself controlled by the "microscopic" 
structural parameters n and $. 

Our computation lays on the identification of the 
elastic energy of the fiber when a force F and a torque 
M are applied at its ends and the sum of the elastic 
energies of its linkers (bend and twist energies) in such 
a constrained configuration. This identity follows from 
the absence of internucleosomal interactions and to our 
focus on the linear response regime (no modification of 
the fiber architecture at this order); it amounts to re- 
strict to the elastic contribution originating from linker 
elasticity within the relaxed chromatin structure. The 
first step is then to relate the "macroscopic" stresses 
(F, M) to the distribution of stresses along linker DNA. 
We denote fj(s) the force and fhj{s) the torque ex- 
erted at the point Pj(s) of linker j (with arclength s) 
by the upstream part of the fiber. We focus on the uni- 
versal, rotationally symmetric behaviour of the fiber. 
Writing standard equilibrium equations of spring me- 
chanics, and using the linear response hypothesis to 
simply sum up the effects of the different stresses ap- 
plied to the fiber (stretching force, bending torque and 
torsional torque), we obtain: 



/iOO = f 

m 3 {s) = M — [Oj(s)Pj(s)} A F 



(1) 



where Oj(s) denotes the orthogonal projection of Pj(s) 
onto the SH axis A. It gives the relation between the 
global stresses F and M exerted at the fiber ends and 
the local stresses fj(s) and rhj(s) experienced at the 
linker level, at each point of the DNA path. The term 
[Oj(s)Pj(s)] A F reflects the involvement of the fiber 
architecture in the expression of the local torque fhj{s). 
This set of equations is the core of our analysis since it 
allows to perform explicitly the change of description 
level, namely passing from the DNA level to the fiber 
level and reciprocally. 

As we restrict to the description of harmonic elas- 
ticity, the coefficients of F and M in the elastic ener- 
gies of the linker will be computed within the relaxed 
SH. We underline that we do not need to compute the 
constrained shape of the linkers to describe the linear 



response of the fiber to applied force and torque. The 
cylindrical symmetry of the fiber ensures that F is di- 
rected along the SH axis A at equilibrium. M decom- 
poses into a bending torque Mb and a torsional torque 
M t according to M = M b + M t A. The well-established 
elastic energy of a linker |J can then be expressed 
as a function of F, M t and Mb thanks to (|l]); it yields 
a quadratic expression 



a s F 2 + a t Mt 



2a t sM t F + a b MZ 



(2) 



which is to be identified with the elastic energy 
£fiber(-F 1 , Mt, Mb) of a length D of chromatin fiber, 
where D is the distance between two successive nu- 
cleosomes measured along the SH axis (projected 
distance). Plugging in the quadratic expression 
£\inkBv{F, Mt, Mt) the linear response ansatz relating 
the stresses F, Mt and Mb experienced by the fiber 
and its strains u (relative extension), f2 (twist rate) 
and p (local curvature): 
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(3) 



naturally embeds the fiber and its elastic behaviour 
within a current continuous model known as the "ex- 
tensible worm-like rope" ( EWLR) in which the 
density of elastic energy writes, as a function of the 
canonical variables u, Q. and p: 



(■EWLR 



k B TAg ksTCQ iu 2 

2 1 2 h — +k B TgQ,u 



(4) 



A is the bend persistence length of the fiber, C its twist 
persistence length, 7 its stretch modulus (dimension 
of a force) and g the twist-stretch coupling constant 
(no dimension). Note the expected vanishing of the 
twist-bend and stretch-bend coupling (for symmetry 
reasons JI(| . We may now express the elastic constants 
of the fiber in terms of the coefficients a b , at, a ts and 
a s , obtained as functions of the DNA elastic constants 
and fiber geometry, by a term-wise identification of 
the coefficients of u 2 , Q. 2 , ufl and p 2 in the following 
identity: 



, k B TAg 2 ksTCQ, ju 4 

tflber = D [ 1 1 \-kBlgilU 



f linker [F{U, tt),M t (u, tt),M b (p)] 



(5) 



from which follows analytical expression of the 
elastic constants A, C, g and 7. We present on Fig- 
ure 3 the results for $ = 50° (high salt concentra- 
tion, with histone HI, i.e. physiological conditions) 
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when the linker length is varied. The curves evidence 
a strong sensitivity of elastic coefficients with respect 
to the fiber structure, here controlled by n. This sensi- 
tivity originates in the sensitivity of the SH pitch: ac- 
cumulation of turns in a spring leads to a highly flexible 
accordion-like behavior, ft is to note that the persis- 
tence length A and C are almost identical, both close to 
the SH pitch and smaller than the corresponding per- 
sistence lengths A and C of DNA. More strinkingly, 
the stretch modulus of the chromatin fiber, around 5 
pN, is to be compared to that of DNA (7dna ~ 1100 
pN): whereas DNA can be seen as non extensible, chro- 
matin is highly extensible; in particular, entropic and 
enthalpic elastic regimes (well-separated in the case of 
DNA) will overlap. We finally underline a key feature: 
the twist-stretch coupling constant, reflecting the chi- 
rality of the chromatin fiber, exhibits steep variations 
and sign changes when the linker length is varied. The 
same twisting of the fiber will either condense or de- 
condense it according to the relaxed structure. 



Experimental implications 

A first issue involving fiber elasticity is to describe 
the linear elastic response of the fiber to global stresses, 
i.e. a force and a torque applied at its ends. This issue 
refers to micromanipulations in which a single chro- 
matin fiber is pulled jl7| (and possibly will be twisted). 
Describing the fiber within the EWLR model, yet thor- 
oughly investigated in the context of DNA by Marko 
|p.6| , allows an interpretation of the force-extension 
curves in terms of geometric and mechanical param- 
eters of the fiber. Experimental results of Cui and 
Bustamante Jl7| are in good agreement with our pre- 
dictions (note that their measure rests on relaxation 
curves, for which there is no interaction between nucle- 
osomes, consistently with our modeling). Indeed, they 
measured a bend persistence length ^4 ~ 30 nm and a 
stretch modulus 7 » 5 pN, which are obtained simul- 
taneously in our modeling for n » 42, or any integer 
multiple of DNA pitch; we there recover the nucleo- 
some phasing observed in biological experiments 

Our study should now be supplemented by investi- 
gating the effect of a dispersion in the values of n and 
<f>. It is likely that the dispersion will induce a kind of 
averaging. Continuity arguments lead us to expect the 
sensitivity of the structural and elastic properties to re- 
main - although possibly smoothed out - together with 
its possible relevance in biological regulation. Work is 
in progress on this question. 

To circumvent this limitation, we suggest to take 
advantage of the interplay between linker intercalation 
and the chromatin fiber mechanics. The twist change 
At achieved in the linker by intercalator molecules (for 



instance with ethidium bromide in vitro) amounts to 
shift the linker length by An = 10.6 x Ar/2vr. Plot- 
ting the various elastic constants as a function of An 
should provide similar curves as those shown on Figure 
3, however smoothed out. We therefore suggest that it 
should be possible to scan experimentally these curves 
by varying continuously the concentration of intercala- 
tors added in the buffer when performing pulling fiber 
experiments. 
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Figure 3: Elastic constants of the chromatin fiber 
versus linker length n (in bp) for $ = 50°. A similar 
sensitivity when n varies is observed for other values 
of $, or when varying $ at fixed n. 
(a) Bend persistence length A (solid red line) in nm 
and stretch modulus in pN (dashed blue line); here is 
shown the effective modulus ~/ e ff = 7— fcaTq jC which 
is involved in force-extension curves |l(| jll) hence di- 
rectly measurable. Our computations predict that A 
and the twist persistence length C are almost equal, 
both close to the pitch V of the SH (Figure 2). Ex- 
perimental measurements of Cui and Bustamante 
yield a value about 30 nm for A and 5 pN for 7 e // 
(horizontal lines). 
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(b) Twist-stretch coupling constant g (no dimension) . 
Note that the vanishing of g corresponds to a change 
of chirality for the fiber. 

Biological relevance 

Another issue, directly relevant to the in vivo func- 
tioning of chromatin, is to investigate both the effect 
of DNA-protein binding on the fiber structure, and the 
effect of fiber stresses on the linker DNA-protein bind- 
ing. Our theoretical approach gives a framework to 
such studies by relating the structure and mechanics 
at the DNA scale and those at fiber scale. 

On one hand, our approach allows to describe the 
response of the fiber to local, internal stresses as those 
created by intercalators, groove-binding proteins, or 
any induced change in the fiber assembly at the linker 
scale. Mechanical sensitivity is likely to provide effi- 
cient switches for processes involving a conformational 
transition of the fiber, for instance the decondensation 
required prior transcription or, at the opposite, fiber 
compaction currently associated with gene silencing 
Jl8[ . A possible tuning mechanism lays on the vari- 
ation of <3>, controlled in particular by the presence 
of linker histone, salt concentration and histone tails 
binding affinities. 

On the other hand, protein-DNA interactions in- 
volved in gene-regulated subcellular events take place 
within chromatin fiber. In consequence, they are likely 
to be influenced, if not directly regulated, by chromatin 
structure and by the stresses generated in the linkers 
when the chromatin fiber is somehow constrained. In 
particular, we expect such an involvement of chromatin 
structure and multilevel mechanical properties to be at 
work in epigenetics and gene imprinting [ju^ . 

Conclusion 

The chromatin scale is precisely the scale of nanome- 
chanics: at this scale, we expect a strong and direct in- 
terplay between the biological functioning, monitored 
by various enzymes, and the mechanical properties of 
the chromatin fiber. Our complete and quantitative 
analysis gives clues to discuss how the elastic prop- 
erties of the chromatin fiber might at the same time 
favour DNA compaction into the chromosomes and al- 
low local and controlled decondensation of chromatin 
involved in gene expression. 
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